Generation of microwave radiation in planar spin-transfer devices 
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Current induced precession states in spin-transfer devices are studied in the case of large easy plane 
anisotropy (present in most experimental setups). It is shown that the effective one-dimensional pla- 
nar description provides a simple qualitative understanding of the emergence and evolution of such 
states. Switching boundaries are found analytically for the collinear device and the spin-flip tran- 
sistor. The latter can generate microwave oscillations at zero external magnetic field without either 
special functional form of spin-transfer torque, or "field-like" terms, if Gilbert constant corresponds 
to the overdamped planar regime. 

PACS numbers: 85.75.-d, 75.40.Gb, 72.25. Ba, 72.25. Mk 
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Spin-polarized currents are able to change the mag- 
netic configuration of nanostructures through the spin- 
transfer effect proposed more than a decade ago [l|, Q. 
Intensive research is currently directed at understanding 
the basic physics of this non-equilibrium interaction and 
designing magnetic nanodevices with all-electric control. 

Initial spin-transfer experiments emphasized the cur- 
rent induced switching between two static configurations 
Q. Presently, the research focus is broadening to in- 
clude the states with continuous magnetization preces- 
sion powered by the energy of the current source @, 0, • 
Spin-transfer devices with precession states (PS) serve as 
nano-generators of microwave oscillations with remark- 
able properties, e.g. current tunable frequency and ex- 
tremely narrow linewidth A particular issue of 
technological importance is the search for systems sup- 
porting PS at zero magnetic field. Here several strate- 
gies are pursued: (i) enginee ring unusual angle depen- 

ii) relying on 
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dence of spin-transfer torque 
the presence of the "field-like" component of the spin 
torque 13 1, (iii) choosing the "magnetic fan" geometry 



14, 15, 



lalrrj. 



PS are more difficult to describe then the fixed equi- 
libria: the amplitude of precession can be large and non- 
linear effects are strong. As a result, information about 
them if often obtained from numeric simulations. Here 
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FIG. 1: Planar spin-transfer devices. Hashed parts of the 
devices are ferromagnetic, white parts are made from a non- 
magnetic metal. 



we study PS in planar devices |18| u sing the effective one- 
dimensional approximation [19 . l2Ci [21 1 which is relevant 
for the majority of experimental setups. It is shown that 
planar approximation provides a very intuitive picture al- 
lowing to predict the emergence of precession and subse- 
quent transformations between different types of PS. We 
show that PS in devices with in-plane spin polarization 
of the current can exist at zero magnetic field without 
the unusual properties (i),(ii) of the spin-transfer torque. 

A conventional spin-transfer device with a fixed polar- 
izer and a free layer (Fig. [I]) is considered. The macrospin 
magnetization of the free layer M = Mn has a constant 
absolute value M and a direction given by a unit vector 
n(t). The LLG equation reads: 
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Here 7 is the gyromagnetic ratio, E(n) is the magnetic 
energy, a is the Gilbert damping constant, s is the spin- 
polarizer unit vector. The spin transfer strength u(n) is 
proportional to the electric current / [H, H3] • In general, 
it is a function of the angle between the polarizer and the 
free layer u(n) = /[(n • s)] /, with the function /[(n • s)] 
being material and device specific 22|, |23j, \24 1 . 

The LLG equation can be written in terms of the polar 
angles (9, <j>) of vector n. Planar devices are characterized 
by the energy form E = (Kx/2) cos 2 9 + E r (9, (j)) with 
K\_ \E r \. The first term provides the dominating easy 
plane anisotropy and ensures that the low energy motion 
happens close to the 9 = ir/2 plane. The residual energy 
E r has an arbitrary form. The smallness of 89 — 9(t) — 
n/ 2 allows to derive a single effective equation on the 
in-plane angle <f>(t) b y p erforming the expansion in small 
parameter \E r \/K± [2l|. For time-independent current 
and polarizer direction s one obtains: 
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where lo±_ = jK±/M. General expressions for a e ff{4>) 
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FIG. 2: (Color online) Evolution of effective energy profile and 
stable solutions with spin-transfer strength (graphs are shifted 
up as u becomes more negative) for a device with collinear 
polarizer. Left: low-field < h < u)u regime. Right: high- 
field h > aii regime. Evolution stage (3) is missing in the high- 
field regime due to the absence of the second energy minimum. 
The red parts of the energy graphs mark the a e f / < regions. 
Filled and empty circle gives represent the effective particle. 



and E e ff((f)) for arbitrary function E r (8,(j)) and polar- 
izer direction s are given in Ref. l2ll In a special case 
frequently found in practice the polarizer s is directed 
in the easy plane at the angle <fi s , and the residual en- 
ergy satisfies (dE r /d8) g=7r /2 — 0; i- e - does not shift the 
energy minima away from the plane. We will also use 
the simplest form /[(n • s)] = const for the spin transfer- 
strength. A more realistic function can be employed if 
needed. With these assumptions 2l| : 
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Equation @ has the form of Newton's equation of mo- 
tion for a particle in external potential E e ff{4>) with a 
variable viscous friction coefficient a e //(0). The advan- 
tage of such a description is that the motion of the effec- 
tive particle can be qualitatively understood by applying 
the usual energy conservation and dissipation arguments. 
In the absence of current, the effective friction is a posi- 
tive constant, so after an initial transient motion the sys- 
tem always ends up in one of the minima of E r (w/2,cj>). 
When current is present, effective friction and energy are 
modified. Such a modification reflects the physical pos- 
sibility of extracting energy from the current source, and 
leads to the emergence of the qualitatively new dynamic 
regime of persistent oscillations. These oscillations of <f> 
correspond to the motion of n along the highly elongated 
(58 <C 1) closed orbits (see examples in Fig. [3J inset), 
i.e. constitute the limiting form of the precession states 



[H H S Ell m spin-transfer systems. 

To illustrate the advantages of the effective parti- 
cle description, consider a specific example of PS in 
the nanopillar experiment where E r is an easy axis 
anisotropy energy with magnetic field H directed along 
that axis, E r (4>) = (K\\/2) sin 2 - HM cosoV The po- 
larizer s is directed along the same axis with <fi s = 
(collinear polarizer). With the definitions wn = r yK a /M, 
h = jH, the effective energy becomes 
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with a)|| = u>ii + u 2 /ujj_. Effective energy profiles are 
shown in Fig. [21 For low fields, \h\ < W||(u), the minima 
at 4> = 0, 7r are separated by maxima at ±(j) m (h). 

According to Eq. ([3]), the effective friction can become 
negative at (f> = or cf> = tt at the critical value of spin- 
transfer strength |m| = u\ = auj±/2. If this value is 
exceeded, the position of the system in the energy min- 
imum becomes unstable. Indeed, the stability of any 
equilibrium in one dimension depends on whether it is 
a minimum or a maximum of E e ff and on the sign of 
a e ff at the equilibrium point. Out of four possible com- 
binations, only an energy minimum with a e ff > is 
stable. A little above the threshold, a e ff is negative in a 
small vicinity of the minimum where the system in now 
characterized by negative dissipation. In this situation 
any small fluctuation away from the equilibrium initiates 
growing oscillations. As the oscillations amplitude ex- 
ceeds the size of the a e / / < region, part of the cycle 
starts to happen with positive dissipation. Eventually 
the amplitude reaches a value at which the energy gain 
during the motion in the a e f / < region is exactly com- 
pensated by the energy loss in the a e f / > region: thus 
a stable cycle solution emerges (Fig. [21 profile (2)). 

The requirement of zero total dissipation means 
that an integral over the oscillation period satisfies 
J a e f f(4>)((f>) 2 dt = 0. In typical collinear systems 25| 
Gilbert damping satisfies a « 0.01 <C W(JuJuj± ~ 0.1 <C 
1, hence the oscillator ©J® operates in the lightly 
damped regime. In zeroth order approximation the fric- 
tion term in ^ can be neglected, and a first integral 
(p 2 /(2uj±) +E e ff = Eq exists. Zero dissipation condition 
can be then approximated by 



OeffWy/Eo-Eeffitfi) d<f> = Q, 



(5) 



with </>i,2 (u) being the turning points of the effective par- 
ticle trajectory, and Eq = E e ff(4>i) — E e ff(4>2)- Since 
the integrand of © is a known function, the formula 
provides an expression for the precession amplitude. 

Consider now the low positive field regime < h < u>\\. 
At u = —ui the parallel configuration becomes unstable 
and a cycle emerges near the 4> = minimum. As u 
is made more negative, the oscillation amplitude grows 
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until eventually it reaches the point of energy maximum 
at u = —U2- Equivalently, the effective particle starting 
at the energy maximum — cf> m is able to reach the other 
maximum at +(f> m (Fig. [21 left, (3)). Above this thresh- 
old the particle inevitably goes over the potential hill and 
falls into the <fi — n minimum which remains stable since 
a eff( 7T ) > holds for negative u. In other words, the cy- 
cle solution with oscillations aroung cf) = ceases to exist. 
At even more negative u the third threshold is reached 
when the effective particle can complete the full rotation 
starting from the energy maximum (Fig. [21 profile (4)). 
Below u = — M3 a new PS with full rotation emerges. In 
the high- field regime h > wn the evolution of the pre- 
cession cycle is similar (Fig. [21 right), but stage (3) is 
missing since there is no second minimum. The thresh- 
old u — —U2 separates the finite oscillations regime and 
the full-rotation regime. 

Thresholds Ui can be obtained analytically from ([5]) by 
substituting the critical turning points 4>\^2 listed above: 
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The corresponding switching diagram is shown in Fig. [3] 
(cf. numerically obtained Fig. 12 in Ref. I2H). It shows 
that different hysteresis patterns are possible depending 
on the trajectory in the parameter space. 

PS in the low field regime was discussed analytically in 
an unpublished work [261 ] . However, since a conventional 
description with two polar angles was used, the calcula- 
tions were much less transparent. Numeric studies of the 
PS were performed in Refs. @,[25| after the experimental 
observation [?[ of the current induced transition between 
two PS in the high field regime. They had shown that 
indeed the low-current precession state PSi has a finite 
amplitude of ^-oscillations, while the high-current state 
PS2 exhibits full rotations of <f> (Fig. [31 inset). 

Next, we consider the cycle solutions in a device called 
a spin-flip transistor [3, [27j . It differs form the setup 
studied above in the polarizer direction, which is now 
perpendicular to the easy axis with <f> s = 7 r/2 . No exter- 
nal magnetic field is applied. In this case 21 1 
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with u>\\ = uju — u 2 /u>j_. As the spin-transfer strength 
grows, the behavior of the system changes qualitatively 
when W|| or a e //|±7r/2 change signs at the thresholds ui — 
±^/W||Cl>_i_ and U2 — ±auj±/2. In accord with previous in- 




FIG. 3: (Color online) Switching diagram of a device with 
collinear polarizer. The u-axis direction is reversed for the 
purpose of comparison with Refs. 0, I25I The parts of the 
diagram not shown can be recovered by a 180-degree rotation 
of the picture. Stable directions in each region are given by 
small arrows, the precession states are marked as PSi,2- The 
large arrow shows the polarizer direction. Inset: schematic 
trajectories of the PSi,2 states on the unit sphere. 



vestigations 2, 28| at |u| > u\ the <f> — 0,7r energy min- 
ima are destabilized and the parallel state <fi = sgn[u]0 s 
becomes stable. Surprisingly, for a > a* = 2y^cJ^JuT[ a 
window u\ < u < U2 of stability of antiparallel configu- 
ration^ = —sgn[u](j) s , opens (Fig. [4]). As discussed in 



Ref. [2l|, the stabilization of the antiparallel state hap- 
pens as the spin-transfer torque is increased in spite of 
the fact that this torque repels the system from that di- 
rection. At u = U2 the antiparallel state turns into a cycle 
(Fig.[H low right panel) which we will study here. Above 
the U2 threshold the amplitude of oscillations grows until 
they reach the energy maximum at u = U3 and the cy- 
cle solution disappears. Although a is not small, U3 can 
still be determined from Eq. ([5]) because a e // is small 
when u is close to U2- Calculating the integral in ([5]) 
with 4>i ,2 = — 7r, we get 



U3 



2 

— aujj_ 



1.27u 2 . (11) 

The usage of approximations |5j) . (fTTj) is legitimate for 
a > 2a» where a e ff(u3) <C tJ&\\(u3)/u)±_ holds. For 
smaller values of a numeric calculations are required. 
They show the existence of a stable cycle down to a = 
0.8 a* where the stabilization of the antiparallel state is 
impossible. For a <C a* and u > u\ the strong negative 
dissipation regime is realized, |a e //| 3> yjui\\/u>±.. Nu- 
meric results show that the amplitude of the oscillations 
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FIG. 4: (Color online) Switching diagram of a spin-flip tran- 
sistor. The u < part of the diagram can be obtained by 
reflection with respect to the horizontal axis. In each region 
stable directions are given by small arrows, precession state 
is marked by PSi. The large arrow shows the polarizer di- 
rection. Threshold 113(a) is sketched as a dashed line where 
approximation (|ll|l is not valid. Lower panels: the evolution 
of effective energy and trajectories (graphs are shifted up with 
growing u) at a « a, (left) and a > a* (right). The red 
part of the energy graph marks the a e ff < region. Effective 
particle is shown by filled and empty circles. 



induced by negative dissipation is so big that the effective 
particle always reaches the energy maximum and drops 
into the stable parallel state (FigdJ low left panel). We 
conclude that the line 113(a) crosses the u = u% line at 
some point and terminates there. 

As for the full-rotation PS, one can show analytically 
that it does not exists in the small dissipation limit at 
a > 2a*. Numerical simulations do not find it in the 
a < a*, u > Ui regime either. 

In conclusion, we have shown that the planar effective 
description can be very useful for studying precession so- 
lutions in the spin transfer systems. It was already used 
to describe the "magnetic fan" device with current spin 
polarization perpendicular to the easy plane 20]. Here 
the switching diagrams were obtained for the spin polar- 
izers directed collinearly and perpendicular to the easy 
direction within the plane. In collinear case we found an- 
alytic formulas for the earlier numeric results, while the 
study of precession solutions in the perpendicular case 
(spin-flip transistor) at large damping is new. The lat- 
ter shows the possibility of generating microwave oscil- 



lations in the absence of external magnetic field with- 
out the need to engineer special angle dependence of the 
spin-transfer torque or "field-like" terms. The inequality 
a > 2* s /u>\\/u)±_ required for the existence of such oscil- 
lations can be satisfied by either reducing the in-plane 
anisotropy, or increasing a due to spin-pumping effect 
(2^ . Most importantly, the effective planar description 
allows for qualitative understanding of the precession cy- 
cles and makes it easy to predict their emergence, sub- 
sequent evolution, and transitions between different pre- 
cession cycle types. E.g., in the systems with one re- 
gion of negative effective dissipation, such as considered 
here, it shows that no more then two precession states, 
one with finite oscillations and another with full rota- 
tions, can exist. Numerical approaches, if needed, are 
then based on a firm qualitative foundation. In addition, 
numerical calculations in one dimension are easier then 
in the conventional description with two polar angles. 
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